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Abstract 

H ! 

Three methods to describe collective motion, Random Phase Approximation (RPA), 
■ Wigner Function Moments (WFM) and the Green's Function (GF) method are compared 

in detail and their physical content analyzed on an example of a simple model, the har- 
monic oscillator with quadrupole-quadrupole residual interaction. It is shown that they 
give identical formulae for eigenfrequencies and transition probabilities of all collective 
excitations of the model, including the scissors mode, which is the subject of our special 
attention. The exact relation between the RPA and WFM variables and the respective 



dynamical equations is established. The transformation of the RPA spectrum into the one 



of WFM is explained. The very close connection of the WFM method with the GF one is 
demonstrated. The normalization factor of the "synthetic" scissors state and its overlap 
with physical states are calculated analytically. The orthogonality of the spurious state 
to all physical states is proved rigorously. A differential equation describing the current 
lines of RPA modes is established and the current lines of the scissors mode analyzed as 
a superposition of rotational and irrotational components. 
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1 Introduction 



The full analysis of the scissors mode in the framework of a solvable model (harmonic oscillator 
with quadrupole-quadrupole residual interaction (HO+QQ)) was given in [1]. Several points 
in the understanding of the nature of this mode were clarified: for example, its coexistence 
with the isovector giant quadrupole resonance (IVGQR), the decisive role of the Fermi surface 
deformation, and several things more. The Wigner Function Moments (WFM) method was 
applied to derive analytical expressions for currents of both coexisting modes, their excitation 
energies, magnetic and electric transition probabilities. Our formulae for energies turned out 
to be identical with those derived by Hamamoto and Nazarewicz [2] in the framework of the 
RPA. However, little details are given in that reference of letter form. In [3] we investigated 
the relation between formulas for transition probabilities derived by the two methods. It was 
shown there, that also these formulas are identical. This coincidence motivated us to undertake 
a systematic comparison of the two approaches and to understand the connection and differences 
between them. This is the goal of the present paper. One of the important subjects of this 
comparison are the current distributions. The WFM method, a priori, can not give the exact 
results, because it deals only with integrals over the whole phase space. It would therefore be 
very interesting to evaluate the accuracy of this approximation by comparing the results with 
the currents obtained from RPA. Unfortunately, even for this simple model (HO+QQ) it is 
impossible to derive in RPA the closed analytical expressions for currents of the scissors mode 
and IVGQR. That is why we consider in addition the Green's Function (GF) method, which 
allows one to find explicit expressions for the currents directly. 

The HO+QQ model is a very convenient ground for this kind of investigation, because all 
results can be obtained analytically. There is no need to describe the merits of the RPA or of the 
GF method - they are very well known [4]. It is necessary, however, to say a few words about 
the WFM. Its idea is based on the virial theorems of Chandrasekhar and Lebovitz [5]. Instead of 
writing the equations of motion for microscopic amplitudes of particle hole excitations (RPA), 
one writes the dynamical equations for various multipole phase space moments of a nucleus. 
This allows one to achieve a more direct physical interpretation of the studied phenomenon 
without going into its detailed microscopic structure. The WFM method was successfully 
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applied to the study of isoscalar and isovector giant multipole resonances and low-lying collective 
modes of rotating and nonrotating nuclei with various realistic forces [6]. The results of WFM 
were always very close to similar results obtained with the help of RPA. In principle, this should 
be expected, because the basis of both methods is the same: Time Dependent Hartree-Fock 
(TDHF) theory in its small amplitude approximation. On the other hand it is evident, that 
they are in general not equivalent, because the truncation scheme is different. The detailed 
analysis of the interplay of the two methods turns out to be useful also from a "practical" point 
of view: first and most importantly it allows one to obtain additional insight into the nature 
of the scissors mode; second we find new exact mathematical results for the considered model. 

The paper is organized as follows. In Section 2 we recall the principal points of the WFM 
formalism and give a summary of the key results of [1] obtained by applying this method 
to the HO+QQ model. The same model is considered in Section 3 in the frame of RPA: the 
formulae for eigenfrequencies, electric and magnetic transition probabilities of the scissors mode 
are derived, the "synthetic" scissors and spurious state are analyzed. The exact interrelation 
between the RPA and WFM methods and between their variables is established in Section 4. 
Section 5 is dedicated to the GF method. The three methods are applied to derive analytical 
formulae for lines of currents in Section 6. The mutual interplay of the three methods is 
discussed in the conclusion. Various mathematical details are given in Appendices A and B. 



2 The WFM method 

The basis of the method is the TDHF equation for the one-body density matrix p r (ri,r 2 ,t) = 
(ri|p r (t)|r 2 > : 



ih 



~dt 



H T , P T 



where H T is the one-body self-consistent mean field Hamiltonian depending implicitly on the 
density matrix and r is an isotopic spin index. It is convenient to modify equation (1) intro- 
ducing the Wigner transform of the density matrix 



/s s 
d 3 s exp(-ip • s/h)p T (r + - , r - - , t) 



(2) 
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and of the Hamiltonian 



H T w (r,p) 
Using (2,3) one arrives [4] at 



d 3 s exp(— ip ■ s/h)(r + 



H 1 



r . 



(3) 



df 



dt 



2 \h r 
— sin < — 
h I 2 



( v)-.(v^-(vr-(v/j|^A (4) 

where the upper index on the bracket stands for the function on which the operator in these 
brackets acts. It is shown in [6, 7], that by integrating equation (4) over the phase space {p,r} 
with the weights x^x^ . . . Xi k pi k+1 . . •Pi n _ 1 Pi„, where k runs from to n, one can obtain a closed 
finite set of dynamical equations for Cartesian tensors of the rank n. Taking linear combinations 
of these equations one is able to represent them through irreducible tensors, which play the role 
of collective variables of the problem. However, it is more convenient to derive the dynamical 
equations directly for irreducible tensors using the technique of tensor products [8]. For this it 
is necessary to rewrite the Wigner function equation (4) in terms of cyclic variables 



df 



i 



(5) 



-1 



with 



+i 



1,0 . d s 



V2 y dx 1 



dx-i 



_d_ 

dx* 



1,0 .8 
V - 1 = T2 { dx~r l 



8X 7 



r+i = --j=(xi + ix 2 ) , 



r 



X3 , r-i = -^( x i -ix 2 ) 



and the analogous definitions for V+ x , Vq , , and p +1 , p , p_ x - The required 

equations are obtained by integrating (5) with different tensor products of r a and p a . Here we 
consider the case n — 2. 



2.1 Model Hamiltonian, Equations of motion 

The microscopic Hamiltonian of the model, harmonic oscillator plus separable quadrupole- 
quadrupole residual interaction is given by 

1 

H=-2 



Z N 



H = E(£ + ^ 2r *) + « E E E 



1 



* 3 
N 



+ 7> K E (- i nE?2-/ 1 N?2 (i (r 3 ) + E?2^(r 1 )g 2M (r 3 )}, 

M=-2 v£j v£j 



(6) 



where the quadrupole operator g 2M = yj T&k / 5 r 2 Y 2 ^ and N, Z are the numbers of neutrons and 
protons, respectively. The mean field potential for protons (or neutrons) is 

V T (r,t) = \muV + £ (-1)^(0(^(0, (7) 

where Z£ = kQ^ + kQ?.^, = + kQ^ an d the quadrupole moments Ql^it) are 

defined as 

QZp(t) = /d{p,r}^(r)r(r )Pl t) 

with / d{p, r} = 2(27r/i)~ 3 / d 3 p / d 3 r, where the factor 2 appears due to summation over spin de- 
grees of freedom. To simplify notation we omit spin indices, because we consider spin saturated 
system without the spin-orbit interaction. 

/ ]_5 

Substituting spherical functions by tensor products r 2 Y 2fJi = \ — r\ , where 

V %ix 

<T,V 

and Ci£ lv is the Clebsch-Gordan coefficient, one has 

V T = \mu^ + Y,(- l Y Z l-A- (8) 

Here 

Rl,(t) = J d{p,v}r^r(v,p,t). (9) 

Integration of equation (5) with the weights r\ , = {r (g) and p| yields the 

following set of equations [1] : 

^-Rla-— L lu = 0, A = 0,2 
dt v m v 

4^-^ + ma; 2 i^-2v^f:v^+l{ffi}(^)^ = 0, A = 0,1, 2 

tit /7J -_Q 

^ + 2m^ M -4v^E v/27 + T{^}(^J)a, = 0, A = 0,2 (10) 

where {^{j is the Wigner 6j-symbol and the following notation is introduced 

Px,(t) = Jd{p,r}p^r(r,p,t), L^{t) = Jd{p,r}(rp) x ,r(v,p,t). 



By definition q 2lI = VQr^, Q T 2yu = VQR^, Rq = -Ql /\/3 with Q T m = N T (r 2 ) being the mean 
square radius of neutrons or protons. The tensor L\ v is connected with angular momentum by 
the relations L{ = ^5/3, L\ ±l = =F HI)- 

We rewrite equations (10) in terms of the isoscalar and isovector variables R Xfl = R 1 ^ + 
R\ni -^A/i = R\ l j,~Rxf i (and so on) with the isoscalar kq = (k+k)/2 and isovector K\ = (k — k)/2 
strength constants. There is no problem to solve these equations numerically. However, we want 
to simplify the situation as much as possible to get the results in analytical form giving us a 
maximum of insight into the nature of the modes. 

1) We consider the problem in small-amplitude approximation. Writing all variables as a 
sum of their equilibrium value plus a small deviation 

R»(t) = RZ + K x »(t), P x ,(t) = Pll + V x ,(t), L x ,(t) = 1% + C x ,(t), 

Rx»(t) = RZ + n x »(t), P x »(t) = P e x l + P x ,(t), L x ,(t) = Lll + C x ,(t), 

we linearize the equations of motion in 1Z XfJ ,, V Xfl , C Xll and TZ XfJi , V Xtl , C Xll . 

2) We study non-rotating nuclei, i.e. nuclei with LH = LH = 0. 

3) Only axially symmetric nuclei with i?2±2 = ^2±i = R%±2 = R-2L1 = are considered. 

4) Finally, we take 

Rio — Roo — 0- (n) 

This means that equilibrium deformation and mean square radius of neutrons are supposed to 
be equal to that of protons. 

Due to the approximation (11) the equations for isoscalar and isovector systems are de- 
coupled. Further, due to the axial symmetry the angular momentum projection is a good 
quantum number. As a result, every set of equations splits into five independent subsets with 
H — 0, ±1, ±2. The detailed derivation of formulae for eigenfrequencies and transition probabil- 
ities together with all necessary explanations are given in [1]. Here we write out only the final 
results required for the comparison with respective results obtained in the framework of RPA. 
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2.2 Isoscalar eigenfrequencies 

The isoscalar subset of equations with // = 1 is 

£21 - 2£ 21 /m = 0, 

£21 - V 2 i/m + [mu 2 + 2k (<520 + 2Qog)] ^21 = 0, 
V21 + 2{muo 2 + KoQ^]£ 2 i - 6K Q e 2 lC u = 0, 

£11 = 0. (12) 

— o 

TflUJ 

Using the self-consistent value of the strength constant k = — — — , (see Appendix A) and the 

4Qoo 

4 

standard definition of the deformation parameter Q20 = Qootj ^ we reduce (12) to 

7^2i - 2£ 21 /m = 0, 
£21 - V 2 i/m = 0, 

V21 + 2mo; 2 [(l + -)C 2 i + SCu] = 0, 

Ai=0. (13) 

Imposing the time evolution via e - *™ for all variables one transforms (13) into a set of algebraic 
equations. The eigenfrequencies are found from its characteristic equation which reads 

Q 2 [Q 2 -2ou 2 (1 + 5/3)] = 0. (14) 

The nontrivial solution of this equation gives the frequency of the /i — 1 branch of the isoscalar 
GQR 

n 2 = n 2 s = 2u 2 (l + 5/3). (15) 

Taking into account the relation (A. 7) we find that this result coincides with that of [9]. The 
trivial solution Q = Q = is characteristic of nonvibrational mode corresponding to the 
obvious integral of motion £ n = const responsible for the rotational degree of freedom. This 
is usually called the 'spurious' or 'Goldstone' mode. It is easy to find from (13) that Cu = 0. 
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2.3 Isovector eigenfrequencies 

The information about the scissors mode is contained in the subset of isovector equations with 

K 21 - 2C 21 /m = 0, 

C 2 i — p2i/m + m uj 2 + kQH + AkxQH TZ 2 i — 0, 
P21 + 2[mu 2 + - Ai = 0, 

+ 3^20^21 = 0. (16) 

Supposing, as usual, the isovector constant K\ to be proportional to the isoscalar one, K\ = an Q , 
taking the self-consistent value for k and the standard definition of 5, we reduce (16) to 

P21 - 2£ 2 i/m = 0, 

5 

C21 - P 2 i/m + moo 2 {\ - a)(l + -)P 2 i = 0, 

o 

r 

P21 + 2mcu 2 [(l + -)£ 21 + 5£ n ] = 0, 

£ n - mcu 2 5(l - a)^2i = 0. (17) 

Imposing the time evolution via e - *™ one transforms (17) into a set of algebraic equations with 
the characteristic equation 

VI 4 - 2fiV(2 - a)(l + 5/3) + 4cu 4 (l - a)5 2 = 0. (18) 

Its solutions are 



Q 2 ± = u 2 {2 - a)(l + 5/3) ± V^(2 - a) 2 (l + 5/3) 2 - 4^ 4 (1 - a)5 2 . (19) 

The high-lying solution fi + gives the frequency Qi V of the /i — 1 branch of the isovector GQR. 
The low-lying solution f2_ gives the frequency Q sc of the scissors mode. 

We adjust a from the fact that the IVGQR is experimentally known to lie practically at 
twice the energy of the isoscalar GQR. In our model the experimental situation is satisfied by 
a = —2. Then 



^ = 4c, 2 [l + ^W(l + ^) 2 -^ 2 , Q 2 sc ^^(l+ 5 --J(l+ 5 -r-^). (20) 



The scissors mode energy in the limit of small deformation is 

E sc ^^?uvS, (21) 

which is quite close to the result of Hilton [10]: E sc m VI + 0.66 hu 5. Taking hu = 45.2/ 'A 1 / 3 
MeV (what corresponds to r = 1.15 fm used in [11]), one obtains 

E sc « 55.4 5A~ 1/3 MeV, 

which practically coincides with the result of Lipparini and Stringari [11] (formula (18): E sc ~ 
56 5A~ l l z MeV) obtained with the help of a microscopic approach based on the evaluation of 
sum rules. Both results are not very far from the experimental [12] value: E sc m 665A _1//3 
MeV. 

It is interesting to study the role of the Fermi surface deformation for the formation of 
IVGQR and the scissors mode. Neglecting the variable p2i(t) in (17) we find that the frequency 
of IVGQR (being determined mainly by the neutron-proton interaction) is changed not very 
much: 

n?„ = 2^(1 -q)(1+ 6/3). 

Comparing this formula (for a = — 2 ) with (20) one sees, that in the limit of small deformation 
one obtains Qf v ~ 6uJq instead of Qf v ~ 8uq. One should recall that also for the Isovector Giant 
Dipole Resonance the distortion of Fermi sphere plays only a minor role. 

It is also easy to see that omitting P 2 i(t) in (17), one obtains zero energy for the scissors 
mode independent of the strength of the residual interaction. Thus, the nuclear elasticity 
discovered by G.F.Bertsch [13] is the single origin for the restoring force of the scissors mode. 
So one can conclude that this mode is in its essence a pure quantum mechanical phenomenon. 
This agrees with the conclusion of the papers [14, 15]: classically (i.e., without Fermi surface 
deformation) the scissors mode is a zero energy mode. 

2.4 Linear response and transition probabilities 

A direct way of calculating the reduced transition probabilities is provided by the theory of the 
linear response of a system to a weak external field 

F(t) =Fexp(-ittt) +F^exp(iQt), 
9 



where F = J2f=i fs is a one-body operator. A convenient form of the response theory is e.g. 
given by Lane [17] (see also section 4). The matrix elements of the operator F obey the relation 



| < U \F\0 > | 2 = h lim (Q - Q u )< ip\F\ip > exp(-iftt), (22) 

where |0 > and \ v > are the stationary wave functions of the ground and unperturbed excited 
states; ip is the perturbed wavefunction of the ground state, Vt u = [E v — E )/h are the normal 
frequencies, the bar means averaging over a time interval much larger than Q being the 
frequency of the external field F(t). 

Electric excitations are described by the operator 



F = *& = T,fM> k = er 2 Y 2fM = [3rl, P = eJ£ (23) 

s=l V 07r 



whose expectation value (in accordance with [1]) is 



< vi^iv >= = \p( r ^ - ^). 

Transition probabilities are given [1] by the following formulae 

B{E2) SC - 2| < SC |F 21 |0 > | - — -Qoo _ fifi • (24) 



m ™ - 2 | <?7; ipP|n>P- e ^ 5 n (l + V3)^-2(^) 2 

5(£2) ra _ 2| < tV |F 21 |0 > | - — -Q o _ ^ • (25) 



B(E2) ts = 2| < fc| j&|0 > | 2 = — ^Qoo[(l + */3)n* - 2(u><5) 2 ]/[^] 3 . (26) 

777. 07T 

These three formulae can be joined into one expression by a simple transformation of the 
denominators. Really, we have from (19) 



±04 -^L) = ±(& + -&.) = ±2^(2 -a)*(l + 6/3)* -4&(l-a)P 

= 2fi|-2cu 2 (2-a)(l + (5/3) = 2fi|-(2-a)(cu 2 + cj 2 ). (27) 

Using these relations in formulae (24) and (25) we obtain the expression for the B(E2) values 
valid for all three excitations 

R(F2) -2i<,i^in>i 2 - e ^ 5 n (i + V3)^ 2 -2(^) 2 

B{E2) V - 2| < HF 21 |0 > | - — ^Qo 0nuK _- 2{2 _ a)il + 6/3) y (28) 
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The isoscalar value (26) is obtained by assuming a — 1. 
Magnetic excitations are described by the operator 

z - eh 
F = *% = YJM, /^-^•[rxVW^Mi,, ^ = W (29) 

Its expectation value was calculated in [1] to be 

Then one finds the following expessions [1] for transition probabilities 

B(Ml) sc = 2| < SC |^|0 > | 2 = -4) (30) 

B(M1)i0 " 21 < wlFnl ° > 1 " —— Qo ° 6 nuni-w sc ) »»■ (31) 

Using relations (27) in formulae (30) and (31), we obtain the expression for the B(M1) values 
valid for both excitations 

FKMU -91 <- 7 ,|pp|n^ I 2 - 1_a ™ 2 n x2 g ~ 2(1 + <V3)(Z> 2 2 

TTIUIq 1 /3 x^/^ 

Taking into account the relation Q^ ~ - I -A) , which is usually [16] used to fix 

ri 2 V 2 / 

the value of the harmonic oscillator frequency ojq , we obtain the following estimate for the 
transition probability of the scissors mode: 

£?(M1)T= 2 | < SC |FP |0 > | 2 = ( 3 / 2 / 1/6 A 4 / 3 ^ / i 2 v = 0.042^^, 

167T 

which practically coincides with the result of [11]: B(M1) |= 0.0A3A^ 3 5 fj, 2 N , obtained with the 
help of the microscopic approach based on the evaluation of the sum rules. 

Concluding this section one should mention, that all magnetic and electric modes of the 
considered model satisfy the energy weighted sum rule [1] , the contribution of the spurious (or 
Goldstone) mode being nonzero. It is interesting to compare the contributions of the scissors 
mode and the spurious mode. The scissors mode (for small 5) yields: 

2E SC \ < 8 c\Fl\0 > | 2 ~ -J-e 2 -Qoo5 2 . (33) 

1287T m 
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The spurious mode yields: 

2hQ Q \ < Q \K\0 >\ 2 = ^e 2 -Q 00 —^—. (34) 

8n m 1 + 0/6 

It is seen that the contribution of the spurious mode is approximately 16 times larger than 
the one of the scissors mode. This is a very significant number demonstrating the importance 
of excluding the spurious state from the theoretical results. Indeed, to describe correctly such 
a subtle phenomenon as the scissors mode, it is compulsory to eliminate the errors from spu- 
rious motion whose value can be an order of magnitude larger than the phenomenon under 
consideration. 

3 Random Phase Approximation (RPA) 

In this section we now want to derive the analogous equations for energies and transition 
probabilities from standard RPA theory. RPA equations in the notation of [4] are 



n,j 

^ ' {v,ij mn X n j -\- \5ij5 mn {^a m Cj) -|- Vi nm j\ Y n j^ TiVtY m i. (35) 

According to the definition of the schematic model in [4], the matrix elements of the residual 
interaction corresponding to the Hamiltonian (6) are written as 

Vmjin — KrT'QimQjn 

with Q im =< i|g2i| m > and K nn = k pp = k, K np = R. This interaction distinguishes between 
protons and neutrons, so we have to introduce the isospin indices r, r' into the set of RPA 
equations (35): 

( e m — € i)Xmi + E K TT'QimQ]nXnj + E K rr' QimQnj^nj ~ ^^mii 
n,j,T' n,j,r> 

E Krr'Q2iQt X n, + ( C ™ " OKi + E ^ = (36) 

n,j,T' n,j,T' 

The solution of these equations is 

Q T * Q T *. 

X "" = m^ K - ™ = ~Jia+~eL (37) 



"•till 
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With 6 mi = C m - 4 aIld RT = Er' «Vr'C T '- 

The constant C T is denned as C T = J2n,j(Qj n ^nj + Qnj^Tj)- Using here the expressions for 
X*j and Y£j given above, one derives the useful relation 

C T = 2S T K T = 2S T ]T k tt ,C t \ (38) 

where the following notation is introduced: 

s T = J2\Q T m f F2 % )2 (39) 

with E = hQ. Let us write out the relation (38) in detail 

C n - 2S n (nC n + RC P ) = 0, 

C p - 2S p {RC n + kC p ) = 0. (40) 

The condition for existence of a nontrivial solution of this set of equations gives the secular 
equation 

(1 - 2S n n)(l - 2S p k) - AS n S p R 2 = 0. (41) 

Making obvious linear combinations of the two equations in (40), we write them in terms of 
isoscalar and isovector constants C = C n + C p , C = C n — C p 

C - 2{S n + S p )k C - 2(S n - S^KiC = 0, 

C - 2(S n - S p )k C - 2(S n + S p )k!C = 0. (42) 

Approximation (11) allows us to decouple the equations for isoscalar and isovector constants. 
Really, in this case S n = S p = 5/2; hence, we obtain two secular equations 

1 - 2Sk = 0, or 1 - Sk = Sk (43) 

in the isoscalar case and 

1 - 2Sk x = 0, or 1-Sk= -Sk (44) 

in the isovector one, the difference of both lies in the strength constants only. Having in mind 
the relation K\ = a/to, we come to the conclusion that it is sufficient to analyze the isovector 
case only - the results for isoscalar one are obtained by assuming a = 1. 
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3.1 Eigenfrequencies 

The detailed expression for the isovector secular equation is 



~mi 



The operator Q has only two types of nonzero matrix elements Q m i in the deformed oscillator 
basis. Matrix elements of the first type couple states of the same major shell. All corresponding 
transition energies are degenerate: e m — = h(u x — uj z ) = e . Matrix elements of the second 
type couple states of the different major shells with AiV = 2. All corresponding transition 
energies are degenerate too: e m — = Ti{uj x + uj z ) = e 2 . Therefore, the secular equation can be 
rewritten as 

(46) 



2«i £ 2 -e§ # 2 -e|' 
The sums <2o = X] |2m«| 2 and Q2 = X! I Qmi 1 2 can be calculated analytically (see 

mi(AN=0) mi(AN=2) 

Appendix B): 

Qo = %e , Q 2 = %e 2 . (47) 
Let us transform the secular equation (46) in polynomial form 

E 4 - E 2 [{e 2 + e 2 ) + 2 Kl (eoQ + e 2 Q 2 )] + [e\e\ + 2/€ ie oe 2 (eoQ2 + e 2 Qo)] = 0. 

Using here the expressions (47) for Q , Q 2 and the self-consistent value of the strength constant 
(A.3), we find 

£4 _ E 2 {1 _ a/2)(e 2 + £ 2) + (1 _ a)e 2 £ 2 = Q> 

or 

VI 4 - Vl 2 (2 - a)uj 2 + + (1 - a)ui = 0, (48) 

with the notation u)\ = uj 2 + w 2 and uj 4 : = (cj 2 — c^ 2 ) 2 . This result coincides with that of [2]. 
By a trivial rearrangement of the terms in (48) one obtains the useful relation 

n 2 (n 2 -uj 2 + ) = (i- a )(n 2 uj 2 + -ui). (49) 

Inserting expressions (A.3) for u 2 , u 2 into (48), we find oo 2 + = 2uo 2 (l + 5/3), cut = A5 2 lu 4 and 
reproduce formula (18) for the isovector case 

fi 4 - 2ftV(2 - a)(l + 5/3) + 4lj 4 (1 - a)5 2 = 0. 
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Taking here a = 



1 we reproduce formula (14) for the isoscalar case 
ft 4 - 2fiV(l + 5/3) = 0. 



3.2 B(E2)-factors 



According to [4] , the transition probability for the one-body operator F = f s is calculated 

3=1 

by means of the formulae 

< 0|F> >= Y.Ul m XZ + fX), < "\F T \0 >= J2(f T mi X^ + fLX^). (50) 

mi mi 

Quadrupole excitations are described by the operator (23) with / 2/i = er 2 Y" 2 ^ = e<2, where 
e = eJj^- The expressions for X^ ni , are given by formulae (37). Combining these results 



we get 



<Q\F» 1 \v>=2~eKlY,\Qm 



fp2 _ ((? \2 

mi v \ c mi) 

The constant Cf, is determined by the normalization condition 



21KISI = ~eCl. 



(51) 



X _ ( yr,u* yr,u' _ yT,w yT,i/\ 

"f," /l^A^mi ^mi 1 mi 1 mi )i 



mi,T 



that gives 



1 



{Cl 



V\2 



\Qr ■ 



P 12 



(C n ) 2 |Q n -I 2 

J I -^mi I 



(& P ) 2 " (^L) 2 ] 2 (C P ) 2 (S») 2 [£ 2 - {e mi ) 2 } 2 
The ratio C n /C p is determined by any of the equations (40): 

_ 1 - 2S p k 2S n R 



(52) 



(53) 



2S P R 1 - 2S n n 

Formula (52) is considerably simplified by the approximation (11), when S p = S n = S/2, e v mi = 
e mn Qmi = Qmi- Applying the second forms of formulae (43,44) it is easy to find that in this 
case C n /C p = ±1. As a result, the final expression for B{E2) value is 



-i 



B{E2) V = 2| < Q\F$y\v > | 2 = 2g 2 \QE u K\Y.\Qr 



{El ~ e 



2 \2 



(54) 



With the help of formulae (47) this expression can be transformed into 



B{E2) V 
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8n moj 2 a 2 E u 
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{El -el) 2 {El -e 2 2 ) 2 

{nw + -uif 
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(55) 



At first sight, this expression has nothing in common with (28). Nevertheless, it can be shown 
that they are identical. To this end, we analyze carefully the denominator of the last expression 
in (55). Summing it with the secular equation (48) (multiplied by oj\), which obviously does 
not change its value, we find after elementary combinations 



Denom = Q^; - 2ft>* + uo 2 + uot + - fi£(2 - a)u 2 + + (1 - a)wi] 

= wlnl[2nl - (2 - a)u 2 + \ - u 4 _[2Ql - (2 - a)uo 2 + ] 
= (Qlul - ui)[2Ql - (2 - a)u 2 + }. 

This result allows us to write the final expression as 



B(E2) V = -^——^Qoq- 



(56) 



(57) 



16tt mu 2 uu tl„ [2Q 2 U -{2-a)uj\X 
which coincides with (28) (we recall that u)\ = 2a> 2 (l + 5/3), out = 4<5 2 u> 4 ). By simple trans- 
formations this formula is reduced to the result of Hamamoto and Nazarewicz [2] (taking into 

5 e 2 h 

account, that they published it without the constant factor — Qm)- 

32n mu Q 



3.3 B(Ml)-factors 

In accordance with formulae (29), (50), (37) the magnetic transition matrix element is given 
by 

(i 11 ) imQim ( f ll) mi Qmi 



<0\F?i\v>=K*Y, 



E - e p • 



E+e». ■ (58) 
As it is shown in Appendix B, the matrix element (/n)i m * s proportional to Q? m (formula 

(B.16). So, expression (58) is reduced to 



-v^-^ 2 ) p E 

"2cV5 V ™ 
eh 



Q p Q p * Q p .Q p * 



J^ivrX^v e mi) e mi(^f + e mi) 

K^—y=(uj 2 x -uj 2 ) p E u J2v r p2 mi i( v xTv 
cV5 mi c mi [E u - (e m J \ 

With the help of approximation (11) and the expressions (47) for Q 0: Q 2 we find 
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2 cVE 
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(59) 



(60) 
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Relation (49) and the self-consistent value of the strength constant K\ = otK were used in the 
last step. For the magnetic transition probability we have 

This relation between B(M1) and B(E2) was also found (up to the factor l/(20c 2 )) by 
Hamamoto and Nazarewicz [2]. Substituting expression (57) for B(E2) into (61) we repro- 
duce (with the help of relation (49)) formula (32). 

3.4 "Synthetic" scissors and spurious state 

The nature of collective excitations calculated with the method of Wigner function moments 
is quite easily revealed analyzing the roles of collective variables describing the phenomenon. 
The solution of this problem in the RPA approach is not so obvious. That is why the nature 
of the low-lying states has often been established by considering overlaps of these states with 
the "pure scissors state" [18, 19] or "synthetic state" [2] produced by the action of the scissors 
operator which antirotates proton versus neutron distributions 

s x =M-\{if)ii-{if)i:) 

on the ground state 

\Syn >= S x \0 > . 

In the considered model the overlap of the "synthetic" state with the real scissors mode (and 
with IVGQR) can be calculated analytically. Let us at first modify the definition of the "syn- 
thetic" state. Due to axial symmetry one can use the I* component instead of I£, or any of their 
linear combinations, for example, the /i — 1 component of the magnetic operator F[ , which is 
much more convenient for us. The terms are introduced to ensure the orthogonality of 

the synthetic scissors to the spurious state \Sp >= (J n + / p )|0 >. However, we do not need 
these terms because the collective states \v > of our model are already orthogonal to \Sp > 
(see below); hence, the overlaps < Syn\u > will be free from any admixtures of \Sp >. So, we 
use the following definitions of the synthetic and spurious states: 

\Syn >= AT 1 ^ - F^IO >, \Sp >= (Ffi + F^)\0 > . 
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Let us demonstrate the orthogonality of the spurious state to all the rest of the states \v >. 
As the first step it is necessary to show that the secular equation (41) has the solution E — 0. 
We need the expression for S T (E — 0) = S T (0). In accordance with (39), we have 



S T (E) 



£2_ £ 2 E 2_ e 2 



, ST{0) = - 



Qo Q2V 

^0 e 2 



The expressions for QJ, Q2 are easily extracted from formulae (B.10), (B.ll): 



Qo ~ Qoo 
m 
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So we find 



(62) 



fi^(0) = --QSo 
m 
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where, in accordance with (B.12), 



{ul - u 2 z r = --(kQ p 2Q + ROM, (u 2 x - "IT = --(«gs, + «<&). 

m m 



(63) 



(64) 



Finally, we get 



2S P (0) = 
2S n (0) = 



Q20 



'—, 1-2S>(0)k = 



20 



^^20 + ^QW ^^20 "I" ^QW 

QW 1 oon/n\.. ^^20 

t^Q/20 + ^^2 



L -j r , l-2,S n (0)ft: 

K V20 + 



/ 2 o ^20 "r ^V20 

It is easy to see that substituting these expressions into (41) we obtain an identity; therefore, 
the secular equation has a zero energy solution. 

For the second step it is necessary to calculate the overlap < Sp\ v >. Summing (59) with 
an analogous expression for neutrons, we get 

eh „ ^ o o,^^ \Ql 



< Sp\u > = 



eh 



Q 



■ T 12 T 

771 \ ^ TfT (. ,2 2\r \ ^ I ^mi I t mi 

-J=E V ^K V (U X -U Z ) ^ .)r( E 2_ e 2 W 

V « r mi V raw V v ^mi/ 



(65) 



Applying the algebraical identity 
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and remembering the definition (39) of S T we rewrite (65) as 

efr 

<S P \v> = __£/C("x-",W-^(0)) 
cv5E v ~ 



eh K p 



In accordance with (38) and (53), 



- coins* - s»(o)) + {ui - uins- - s-m^ 



(66) 



Kl 1 - 2S p k 

= -2~s^r- (67) 

Noting now (see formula (63)) that {uo 2 x — loI) t S t (0) = —-^Qlo and taking into account relations 
(64), we find 

<Sp\u> = /? {[(/cQg + RQl)2S* - Ql\ + [(«QS + RQl)2S» - Qg] 1 "J^" } 

{i _ 9 9 p k 
[(2S«Vc - i)gp + 25**05] + [(25V - i)g° + 2^/5^)] 2 ^ n _ 

1 -25V 



/3 \2S P RQ* + (2S n n - l)Ql 



I ^ z K 2S n R 
= p^^{2S n R2S p R- (1-2S u k)(1-2S p k)} = 0, (68) 

3 eh K p 

where (5 = 7= — - and Q2 = g2o- The expression in the last curly brackets coincides 

m cV5 /'-',/ 

obviously with the secular equation (41) that proves the orthogonality of the spurious state to 
all physical states of the considered model. So we can conclude that strictly speaking this is not 
a spurious state, but one of the exact eigenstates of the model corresponding to the integral of 
motion J 11 + I p . The same conclusion was made by N. Lo Iudice [20] who solved this problem 
approximately with the help of several assumptions (the small deformation limit, for example). 

The problem of the "spurious" state being solved, the calculation of overlaps < Syn\u > 
becomes trivial. Really, we have shown that < Oj-F^ + Ff^u >= 0, hence < OjF^lz/ >= 
- < 0\F^\u >. Then < Syn\v >= AT 1 < 0|Ffi - FnW >= 2 ^ < °\ F nW > and 

U 2 = \< Syn\v > | 2 = 2N- 2 B(M1) U . (69) 

The nontrivial part of the problem is the calculation of the normalization factor A/". It is 
important not to forget about the time dependence of the synthetic state which should be 
determined by the external field: 

\Syn{t) >= A/'- 1 [(A P i " F^)e~ m + (F£ - F^e m ]\U > . 
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Then we have 



Af 2 = 2<0|(A P i-F 1 n 1 )t(Ff 1 -F 1 n 1 )|0> 

= 2]T < OKFfi - F?J\ph >< phftFZ - F[\)\0 >= 2]T | < P h\{F^ - F^O > | 2 

= 2^|<^|F 1 7 - 1 |0>| 2 = 2 5:|(/ 1 U,| 2 . (70) 

r,ph T,ph 

With the help of relation (B.16) we find 

A/" 2 = -(— ) 2 ]T (V 1 <P h \ r2Y ^\° > I'V 



Expressions for Qq, <2£, w J? w I are given by formulae (62), (B.12). To get a definite number, 
it is necessary to make some assumption concerning the relation between neutron and proton 
equilibrium characteristics. As usual, we apply the approximation (11), i.e., suppose Q^ = 
Qq , Q20 — Qw- It 1S eas y to check that in this case formulae for uj t x z are reduced to the ones 
for the isoscalar case, namely (A. 3), and <2q = Qo/2, Ql = Q2/2, where Qo and Q2 are given 
by (47). So we get 

8n 2c moj 1 \e e 2 / 2n n 

The estimation of the overlap for 156 Gd with 5 = 0.27 gives N 2 = 34.72/i 2 ^ and U 2 = 0.53 (see 
eq. (69)), that is two times larger than the result of [18] obtained in QRPA calculations with 
the Skyrme forces. The disagreement can naturally be attributed to the difference in forces 
and especially to the lack of pair correlations in our approach. 

In the small deformation limit U 2 = \\J\ ~ 0.6. This is the maximum possible overlap of 
the "pure" (or "synthetic") scissors with the real scissors. Increasing 5 and /or taking into 
account pairing correlations decreases its value, that is confirmed by numerous microscopic 
calculations with various forces [21]. Such small overlap leads inevitably to the conclusion, that 
the original model of counter rotating rigid rotors [22] has not very much in common with the 
real scissors mode, the correct description of which requires the proper treatment of the Fermi 
surface deformation and the coupling with IVGQR. 
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3.5 Superdeformation 

A certain drawback of our approach is that, so far, we have not included the superfluidity 
into our description. On the other hand, our formulae (20, 32) can be successfully used for 
the description of sup erde formed nuclei where the pairing is very weak [2, 22]. For example, 
applying them to the superdeformed nucleus 152 Dy (5 ~ 0.6, Hljo = 41/A 1 / 3 MeV), we get 

E iv = 20.8 MeV, B(Ml) iv = 15.9 /j? N 

for the isovector GQR and 

E sc = A.7MeV, B(M1) SC = 20.0/4 

for the scissors mode. There are not many results of other calculations to compare with. As a 
matter of fact, there are only two papers considering this problem. 
The phenomenological TRM model [22] predicts 

£™~26MeV, B(Ml) iv ~26n 2 N , £ sc ~6.1MeV, B{Ml) sc ~ 22(i 2 N . 

The only existing microscopic calculation [2] in the framework of QRPA with separable forces 
gives 

£ iw ~28MeV, B(Ml) iv ~ 37 ij? n , E sc ~ 5 — 6 MeV, B(M1) 1+ ~ 23/^. 

Here B(M1) 1 + denotes the total Ml orbital strength carried by the calculated K n = 1 + QRPA 
excitations modes in the energy region below 20 MeV. 

It is easy to see that in the case of IVGQR one can speak, at least, about qualitative agree- 
ment. Our results for E sc and B(M1) SC are in good agreement with that of phenomenological 
model and with E sc and B(Ml)i+ of Hamamoto and Nazarewicz. 

It is possible to extract from the histogram of [2] the value of the overlap of calculated 
low-lying 1 + excitations with the synthetic scissors state: | < Syn\l + > \ 2 0.4. The result of 
our calculation U 2 = 0.43 agrees with it very well. So, the comparison of our calculations with 
that of QRPA shows, that we have excellent agreement in superdeformed nuclei and rather 
large disagreement in moderately deformed nuclei. On the other hand it is known [2], that 
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pairing is very weak at the super deformation and becomes important at moderate deformations. 
Therefore, as a consequence the correct treatment of pair correlations is important for an 
accurate description of the scissors mode. This shall be the subject of future work. 



4 WFM versus RPA 

In this section we want to establish the precise relation between the WFM and RPA methods. 
Though it can be guessed that working in both cases with the full configuration space, i.e. 
with a complete set of particle and hole states in the case of RPA and with all moments of all 
ranks in the case of WFM, identical results shall be obtained, it is still important to study the 
detailes of this relation because the behaviour of WFM and RPA under truncation of the space 
turns out to be radically different. The exact relations between the RPA amplitudes X ph and 
Y ph and the respective WFM variables can be established with the help of the linear response 
theory. 

Let us consider the response of the system to a weak external time-dependent field 

W{t) = W exp(-int) + W ] exp(iftt) (73) 
with W = J2k q w kqol.a q . The linear response is then given by [4]: 

k'q' 

where 

Rkq,k'q'{&) = E 



< 0|ajafe|z/ >< i/|ajj/a 9 /|0 > < 0|a|,/a g /|z/ >< u\a q ak\0 > 



is the RPA response function [4] and the index pairs kq and k'q' are restricted to particle hole 
and hole particle pairs. Indices p, q include spin and isospin quantum numbers a and r. For 
the change of the average value of an arbitrary operator we have: 

5<^\F\^>=J2fqkP { k 1 q ) - (74) 

kq 

We now are ready to analyze the WFM variables. The first one is 

R^(t) = 2(2nh)- 3 Jd 3 p |rf 3 rr^r(r,p,t). 
22 



Using here definitions of the Wigner function, the 5-function and the density matrix [4] 

p(var, rW) = £ ^(rWr'^rar) < *|4a,|* > (75) 

kq 

we find 

R T ^(t) = ■^^/rf 3 rr^|rf 3 s |rf 3 pexp(-2p-s/% T (r + |,r-|,t) 
= 2 |d 3 r r^p r (r, r, *) = £ |rf 3 r r^p(r<rr, rar, i) 
= EE /d 3 rr^(rar)0,(rar) < *|4a,|* > 

= E( r iX < *l4a«l* >=< *l E( r A /i )^4« 9 l* > 

fcg fcq 

= <*lE^( a )|tf >=<¥|i^|tf >, (76) 

s=l 

i.e. this is just the ground state expectation value of the operator = E( r s)-V»- In accor- 

8 = 1 

dance with (74) the variation of this variable is 

kq 

= E(< °\K\» > c -- < > c») e ~ int 

V 

+ E(< > <£- < OlAyi/ > c:)e^, (77) 

V 

where 

< i/|W"|0 > _ „ < i/|4a g |0 > _ < 0|W"|i/ > _ < Q\a\a q \v > 

Cv = h(Q -n v ) = \; h(Q - n„) Wkq > Cv = h(Q + n,) = ^ h(Q + Q u ) Wkq - (78) 

Equation (77) demonstrates in an obvious way the structure of the variables 5-Ra^- They 
are linear combinations of the transition matrix elements < OI-Ra^I^ > which are, in turn, linear 



A 

E 

s=l 



combinations of the RPA amplitudes X kq , Y kq . Introducing the notation L Am = '^2(r s p s )\^ and 

A 

= E(PsW we can proceed in a similar way with 

L^it) = 2(27r/r 3 / d 3 p J AM v f(r, p, t) 



8 = 1 
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and 

P A yO = 2(27rhy 3 Jd 3 p jdhplj^v^t). 

Inserting the expressions for SRx^ SL\^ into the first equation of the set (10) we find 

2 - 
-ift^(< 0|i2 A/ > > c„- < ^|i?A M |0 > c„) = — XX< °l L A M k > c^- < v\L Xtl \0 > c„). 

v m V 

It is sufficient to consider only the part with the e~ time dependence. Multiplying this 
equation by (Q — Q u ) and taking the limit Q — > Q u we find the equation 

2 

< OlRxpW >= — < 0\L XtM \u >, (79) 
m 

which can be called as the dynamical equation for the transition matrix element < 0|i?A^|^ >• 
Naturally, in the same way the dynamical equations for transition matrix elements < OILa^I^ > 
and < 0|Pa m |^ > can be extracted from the second and third (linearized) equations of (10). 

Now we can show, that exactly the same dynamical equations for Transition Matrix Ele- 
ments (TME) can be derived from RPA equations. To this end we combine the RPA equations 
(36) with the definition (50) of matrix elements: 

hVlu ^^X fim-^-mi fmiXmi ) ^ ] €mi(j ~im-^ mi fmi^mi ) ^ X fimQjm fmiQmi)' (80) 

mi mi mi 

Taking into account the relations 

^mifim [ft -^o]im> ^mifmi [ft -^o]mi> 

one rewrites this equation as 

«l„ < 0|F> >= ]T{[/\ Hl\ lm XZ + \f, Hl\ ml Y^ + KiflnQZ - fmiQmi)}- (81) 

mi 

The Hamiltonian of the axially deformed harmonic oscillator corresponding to the mean field 
(8) is 

#o T (r) = £{J| + + ^ (eg)r 2 2 ( S )}. (82) 

Let us consider the operator / = \[§r\ x = q 2 \ = Q. Calculating the commutator 

2 

[r%i,Ho\ = ih—(rp) 2 i 
m 
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we find from (81) the following equation 

mi v < o| E(^ 2 i)Ik > = in- E{(M 2 i)L^r + ((rphiY^n 

s=l m mi 

+K tS T(Q t Q t * - Q T Q T * ). (83) 

mi 

Taking into account relations (Q*)i m = (Q) m i and |Q m i| 2 = |Qj m | 2 we find, that the last sum 
in (83) is equal to zero. Applying again formula (50) we write (83) as 

-iSl v < 0\R^\u >= — < 0\L T 21 \i> > (84) 

reproducing equation (79). The dynamical equations for < OlL^I^ > and < Ol-P^lz/ > are 
obtained in a similar way by considering the operators / = (rp)| x and / = respectively. 
As it could be expected, the resulting equations reproduce the corresponding TME equations 
derived from (10) with the help of the limit procedure. 

So, there exists one-to-one correspondence between the set of dynamical equations for WFM 
variables and the set of dynamical equations for Transition Matrix Elements (TME). This 
correspondence makes obvious the fact that both sets have the same eigenvalues. On the other 
hand the TME equations are just linear combinations of the RPA equations. Therefore we 
can conclude that RPA and WFM approaches generate identical eigenvalues. In this sense 
both approaches are equivalent in all aspects. This concerns for instance also the transition 
probabilities. However, for this equivalence to be exact, one needs to work in the full space in 
both approaches, that is in the complete particle hole space in RPA and taking all phase space 
moments of all powers in WFM, a task which can hardly be tackled in general. 

The difference of the two approaches then shows up if truncations of the dimension of the 
equations have to be operated. In RPA one usually solves the equations with a restricted 
number of discrete particle hole pairs, i.e. the dimension of the RPA matrix is finite (in some 
works the RPA equations for finite nuclei are, however, solved in full space, including continuum 
states [23, 24]). The result of such a diagonalisation usually yields a huge number of discrete 
eigenvalues approximating more or less the spectrum one would obtain from a solution in the 
full space. For instance resonances in the continuum (e.g. giant resonances) will be mocked 
up by a bunch of discrete states whose envelope may simulate the full solution. Reducing the 
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dimension of the particle hole space too much may lead to a situation where the full solution 
is only approximated rather badly and in an uncontrolled manner. 

In the WFM method the dynamical equations for Cartesian tensors of the rank n = 2 are 
coupled (by the interaction terms in (5)) with dynamical equations for tensors of the rank 
n — 3, these equations being coupled with the ones for tensors of the rank n = 4 and so 
on up to n = oo. Here one hopes that the essential part of physics is described by a small 
number of the lowest ranks tensors. The hope is based on the assumption that the higher 
rank tensors (moments) are responsible for the more refined details and that neglecting them 
does not appreciably influence the description of the more global physics which is described 
with the lower ranks tensors. This assumption is substantiated in past applications of the 
WFM method to realistic situations with Skyrme forces for the description of collective nuclear 
modes [6, 25, 26]. In those works it has indeed been demonstrated that even with a very limited 
number of low rank phase space moments one can faithfully reproduce the centroid position of 
the collective states. In this sense the WFM method is rather similar to the sum rule approach 
which works, however, only in the cases when practically all strength is exhausted by one state, 
whereas WFM method works also in situations when the strength is distributed among a few 
excitations. From these studies it is then permitted to assume that the inclusion of higher and 
higher rank moments will just give raise to a refinement of the gross structure obtained with the 
low rank tensors. A formal convergence study of this type has been performed in the infinite 
matter case [27] where it was indeed shown that the moment method allows to approach the 
full solution in an optimized way. 

The net result is, that WFM and RPA approximate the exact infinite spectrum into a finite 
number of eigenfrequencies with, however, different convergence. 

An analogous situation occurs with transition probabilities. Let us analyze, for example, 
the expression (77) for the WFM variable SR^(t) = U T Xll (t). Using the definition (78) of c v 
with the external field operator W = R T X „, we find 

xrr (f] f ( < o\K\» >< »\K\° > < H^|o><QK>> \ _«* 

*° (\ < 0|i^> > | 2 | < 0\R T x l\u > \ 2 \ _ iQt 
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The summation limit iV c depends on the method of calculation. In the case of the exact solution 
N c = oo, for RPA N c is usually of the order of several hundredsorthousands, for WFM N c 
usually is not more than around a dozen. Naturally, the eigenvalues VL V and eigenstates \v > are 
different in each case. So, the strength, that in RPA was distributed over hundreds or thousands 
levels, in WFM is concentrated only on several levels, i.e. averaging of levels is accompanied by 
the redistribution of the strength. The variable 5R{ (t) is the quantum mechanical observable, 
so its value should not depend on the basis \v >. Hence, the right hand sides of (85), calculated 
by two methods, should coincide if both methods are mutually consistent. This statement can 
e.g. be checked with the help of sum rules. Generally in RPA sum rules are well fulfilled 
for a sufficiently large particle hole space which in realistic cases can become quite significant, 
whereas in WFM sum rules are generally already well fulfilled even with a small number of low 
rank moments (see e.g. [6, 25, 26, 28]). 

The essential difference between WFM and RPA methods lies in their practical use. The 
RPA equations (35) are constructed in such a way that the increase in dimension does not cause 
any formal problems and finally it is only a question of computer power what dimension can be 
handled. Quite on the contrary the increase of dimension in WFM is a nontrivial task. Beyond 
a certain order of the moments even the reduction of Cartesian tensors to the irreducible 
ones becomes a very difficult task. However, the spirit of WFM is rather to reproduce the 
gross structure of a couple of prominent collective states, a situation which it can handle very 
efficiently However, in cases where there is strong fragmentation, direct diagonalisation of the 
standard RPA equations is more efficient. 

In conclusion WFM and RPA are equivalent when the full particle hole configuration space 
in RPA and the infinite number of moments in WFM are considered. However, under truncation 
of the spaces both methods have different convergence properties. In the general case for WFM 
only a few moments are sufficient to get the correct gross structure of the collective part of the 
spectrum, whereas in RPA one in general must take into account a quite large configuration 
space to produce reasonable results. 
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5 Green's function method 



One of the important subjects of comparing RPA and WFM methods are the current distribu- 
tions. The WFM method, a priori, can not give the exact results, because it deals only with 
integrals over the whole phase space. It would therefore be very interesting to evaluate the 
accuracy of this approximation by comparing it with the exact result. Unfortunately, even for 
the simple model HO+QQ it is impossible to derive in RPA closed analytical expressions for 
currents of the scissors mode and IVGQR. That is why we consider in this section the Green's 
Function (GF) method, which allows one to find explicit expressions for the currents directly. 

Following the paper of H. Kohl, P. Schuck and S. Stringari [29] we will consider at first 
the isoscalar case. Conserving on the right hand side of equation (4) only the first term of the 
sin-function expansion leads to the Vlasov equation 

?l = VH W ■ V p f - V P H W ■ Vf. (86) 

In our case the Wigner transform Hw coincides with the classical Hamiltonian H c . Having in 
mind small amplitude vibrations we have to linearize (86): f — fo + fi, H c — H + H 1 , with 
f being the solution of the time independent equation. The linearized version of (86) is 

^ + V P H ■ VA - Vif • y p fi = S(r, p, t), (87) 

where S(r, p, t) = Vi?i ■ V p / - This equation will be solved with the Green's function method. 
We have 

(4- + V P H • V - Vif • VnG^frp, r'p') = S(r - r'Wp - p')6(t - t') (88) 
ot 

with [29] 

GC-'Vp, r'p') = a[r c (r, p, t'-t)- r']5[p c (r, p, t'-t)- p']6(t - t'), 

where r c (r, p, t' — t), p c (r,p,t' — t) are solutions of classical equations of motion with initial 
conditions r, p. The solution of (87) can be written as 

oo 

/i(r,p,f) = tf+ J dt' j Sr'dyG^ivvy^Siv'^'X) 

— oo 

t 

= /f+ J dt'S(r c ,p c ,t'), (89) 

— oo 
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where f± is the solution of the homogeneous equation. It is obvious that any function of 
variables r c and p c satisfies the homogeneous equation, however it does not play any role at 
resonance and we therefore will omit it in the forthcoming. 

We consider the axially deformed harmonic oscillator with the quadrupole-quadrupole resid- 
ual interaction V res . Therefore the single-particle Hamiltonian is 

tfo = fUy[^ 2 + y 2 )+- 2 V] 

and Hi = V res . 

We are interested in the part of the residual interaction with = 1. In accordance with 
formula (A.l) it can be written as 

V res = -«o[<52i(0<72-i(r) + Q 2 -i(t)q 2 i(r)] = 12K Q Qi(t)[xz + yz\ 

with 

Qi(f) = 2 Jd{p,r}f(r,p,t)xz = 2 Jd{p, r}/(r, p, t)yz 
= 2 /d{p,r}[/ (r,p) + fi{r,p,t)]xz = 2 Jd{p,r}fi(r,p,t)xz. 
With the help of the Thomas-Fermi approximation for the static distribution function 

/ = 9(e F - H ) 

the right hand side of (87) is found to be 

S(r,p,t) = -12— Qi(t)6(e F - H )\p x z + p g x + p y z + p x y]. 
m 

The classical trajectories are determined by the solution of Hamilton equations 
f c ,i = Pc,i — with i — x,y, z. In our case they are 

r ci(t) = Ti cos u{t H — sinwjt, p C i(t) = p% cos ourf — mujiTi sin u{t. 

muji 

Formula (89) then gives 

t 

A(r,p,t) = -Q^5(e F -H ) [ dt'Qi(t') x 

— oo 

x{ — [p x + p y )z [uj + COSUJ + (t' — t)+UJ_ cosc<j_(t / — t)] 
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+ — p z (x + y) [uj + cos iv+(t' — t) — uj- cos cj_ (t' — t) ] 

H (p x + p y )p z [uj+ sinu + (t' — t) — uj_ sinu;_(t' — t)] 

muj x uj z 

—m(x + y)z[u + sinc<j + (t / — t) + uo_ sinc<j_(t / — t)]}, (90) 

where u± = uj x ± uj z . 

So, we have derived a complicated integral equation for the perturbed distribution function 
which may not easily be solved in general. As a matter of fact the analytic possibilities of the 
Green's function method are, without further consideration, exhausted at this point. However, 
expressions (88) and (89) point to the possibility to use the so-called pseudo particle method 
[30], in case the classical trajectories are not known analytically. 

In order to proceed to the evaluation of the eigenfrequencies and transition probabilities 
we again apply the method of moments. Integrating (90) over the whole phase space with 
the weights xz, p x p z , zp x + xp z and zp x — xp z we obtain the following set of coupled integral 
equations 

t 

Qi(f) = J dt'Q 1 {t')[uj + sin u + {t' -t) +w_sinw_(t' -t)}, 

— oo 

t 

Pi(t) = —{3m 2 uj x uj z J dt'Qi(t')[u + smu+(l/ — t) — u>- sinc<j_(t / — i)}, 

— oo 

t 

L^t) = -13m J dt' Qi{t')[uj 2 + cosuj + (t' — t) + uj 2 _ cosuj-(t' — t)], 

— oo 

t 

I y {t) = -pmcu + LU_ J rft'Q 1 (t , )[coscj + (t' - t) + cosw_(^ - *)], (91) 



where 



_ 2n^e F 4 



m 2 ui x ul (27rh) 3 



= 12k J d{p, r}x 2 z 2 5(e F - H ) = J d{p, r} P 2 x p 2 z S(e F - H ) 

and the following notation is introduced 

Pi(*) = 2 Jd{p, r}/i(r, p, t)p x p z , L^t) = 2 Jd{p, r}/i(r, p, t)(zp x + xp z 
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I y (t) = 2 J d{p,r}/i(r,p,t)(zp.,. -xp z ). 



By simple means these equations are reduced to a set of differential equations. At first we 
perform time derivatives of all equations in (91): 

t 

Q x {t) = -13 J dt'Qxit'^l cos u + {t' -t) + cu 2 _ cosct>_(t' - t)}, 



— CO 

t 



P 1 (f) = (5m 2 uj x uj z J dt'Qi (*') [lu 2 + cos lu + (f -t)-uj 2 _ cos lu_ (t' - t)], 

— oo 

Li(t) = -j3m j(c4 +ui)Qi(t) + y rft / gi(t / )[^+sincj + (t' -t) +ulsmu-(t' -t)]J , 

j y (t) = -pmu + u_ |2Qi(t) + y dt'Q 1 (t')[uj + sincj + (t' -t)+U- smu^(t' - . (92) 

Solving (91) with respect of four time integrals (containing sinu;± and cosu;±) we can substitute 
found expressions into (92). We obtain 

Qi(t) = -Mt), 

m 

U{t) = -m(2p+l)(u; 2 x + u; 2 z )Q 1 (t) + -P 1 (t), 

m 

i y (t) = -m(2(3 + l)w+w_Qi(t). (93) 

Due to conservation of the angular momentum the right hand side of the last equation must 
be equal to zero. So we have the requirement 



m 2 ^ x u 3 z (2nh f 
4vr 3 4 4 



2/5 + 1 = 0, or . = - T:7 " 7 • (94) 



With the help of the relation 



A(0 = " = ^Mffl) ^ (95) 



and formulae (A. 3) for a; x , u z the expression for k is reduced to 



m{u 2 x + 2u 2 ) muj 2 
12A(r 2 ) = ~4A{r 2 } 



ntyuj x -r ^ z mm 
K = 10 4/„2\ = "TTT^v ( 96 ) 



31 



which is just the familiar expression for the self-consistent value of the strength constant (see 
Appendix A). This is a rather interesting result, because the well known formula is obtained 
without the usual self consistency requirement [31]. As it is known, in the absence of exter- 
nal fields the angular momentum of any system is conserved. The short range interparticle 
interactions depending on the module of the interparticle distance |rj — rj\ create the scalar, 
i.e. rotational invariant, mean field, which exactly repeats the shape of the nucleus. When 
we imitate the mean field by a rotational invariant function, the angular momentum will be 
conserved independently of the shape of this function due to a pure mathematical reason: an- 
gular momentum operator commutes with a scalar field. If we use the non rotational invariant 
function (as in our case), mathematics does not help and the shape of the function becomes 
important. If the function does not follow exactly the shape of the system, the latter will react 
on this inconsistency as on the external field, that leads to the nonconservation of an angular 
momentum. Therefore the requirement of the angular momentum conservation in this case 
becomes equivalent to the requirement of the self consistency. This is seen very well in the 
method of moments. Integrating equation (87) over the phase space with the weight zp x — xp z 
we obtain the dynamical equation for I y 



Obviously it is the requirement of the consistency between the shapes of the potential and 
the nucleus. In principle this relation is less restrictive than the standard self consistency 
requirement [31]. However, the latter satisfies equation (98) what can be easily checked with 
the help of Appendix A. 

So, finally the set of equations (93) is reduced to 




(97) 



The requirement of the angular momentum conservation gives the following relation 




(98) 



Qi(t) 



l 



m 



Li(t) 



2 



m 




Iy{t) 



0. 



(99) 
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Taking into account the relations between the definitions of variables in (13) and (99) 



Qi = -Re7Z 21 , Pi = -ReV 2 i, L x = -2ReC 2 i 



(which follow from formulae r\ 



21 — 



z(x + iy) and (rp) 



21 - -o 



\[zp x + xp z + i(zp y + yp g )]) and 



I y = 2ReCu, it is easy to see, that the last set of equations is identical to (13). 

With the help of relations (91) the Wigner function (89) can be written in terms of the 
Wigner function moments. Taking into account equations of motion (99) and the time depen- 
dence of variables via e~ %nt (which leads to the equality I y = 0) one finds 



In the case of 5 = it reproduces the result of [29]. 

Having the Wigner function one can calculate transition probabilities in the same way as 
in WFM method. 

Let us consider now the problem with two sorts of particles: neutrons and protons. All vari- 
ables and parameters acquire isotopic index r. The part of the residual interaction with = 1, 
in accordance with formula (8) becomes V{ = Z\{t)\xz + yz] with Z"(t) = 12(/cQi + kQi), 
Zf(t) = 12 (kQ * + kQD and Q\ = fd{p, r}f{(r, p, t)xz. The expression for the Wigner func- 
tion is obtained from formula (90) by changing the factor QKQi(t') by ^Z{(t'). The dynamical 
equations for isovector variables Q l = Q\- Q\, P 1 = P* - Pf, L X = L\- L\, and I y = If - 1? 
can be derived (in approximation (11)) exactly in the same way as the equations for isoscalar 
ones. As it is expected, they coincide with (17). 

As we see, in the considered simple model all results of WFM method are identical to 
that of Green's Function (GF) method. Having in mind also that both methods generate the 
same set of dynamical equations for collective variables (Wigner function moments), one could 
suspect their identity. In general, this is not quite true. The principal difference between the 
two methods is more or less obvious. In the GF method one finds first the formal solution of 
eq. (5) and only afterwards one takes the phase space moments of the found Wigner function 
to obtain the final solution of the physical problem. In the WFM method one takes from the 




3 k 

f3m 2 
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beginning the phase space moments of equation (5) without any attempts to find the "natural" 
expression for the Wigner function. 

The reason of coincidence of all results is quite simple. For the harmonic oscillator with 
multipole-multipole residual interaction of arbitrary rank (multipolarity) the equations of both 
methods can be derived without any approximations - the interaction of the multipolarity n 
generates the set of dynamical equations for tensors (moments) of the rank n. For the GF 
method this is easily seen from formula (89). In the case of the WFM method it is seen very 
well from the structure of equation (87). When one takes the moments of rank n, neither the 
left hand side no the right hand side of this equation can generate moments of rank higher than 
n. The coincidence of results in the case of n = 3 was demonstrated in [32]. 

The power and simplicity of the GF method are restricted by the potentials for which 
the analytical solutions for classical trajectories are known. In the case of realistic forces the 
GF method loses its simplicity and transparency, however the pseudo particle method [30] 
can still be applied. The WFM method does not meet any difficulties and continues to be a 
convenient and powerful tool for the description of the collective motion what was demonstrated 
by calculations with Skyrme forces [6] . For an illustration of this property of the WFM method, 
currents are a good example, because the procedure of their construction with WFM is general 
enough to be used for any type of force (see section 6.1 below and [1]). 



We are interested in the trajectories of infinitesimal displacements of neutrons and protons 
during their vibrational motion, i.e. in the lines of currents. The infinitesimal displacements 
are determined by the magnitudes and directions of the nucleon velocities u(r, t), given by 



6 Flows 




2^{(V-V>(r,r',0} 



V E(( V " V>(nxr, r'ar, t)} r=r , 



^EEW^V^rar) 



? (rcrr)V0;(r ( 7r)} < *|aja 9 |* > 



pq <t,t 



34 



pq pq 

= m<^\j(r)\y>. (100) 
The current density operator J(r) has the standard quantum mechanical definition [4]: 

A ■* A 

i(r) = E j.(r) = ~ E^( r " r s )V s + V s 5(r - r.)] = Y.UM a v 

s=l Zm s=l 

jft _ 
J Pq (r) = -— < p\[5(r - r)V + V5(r - r)]\q > 

= ^EiU^V^var) - ^{var^Uvar)} = 4^-[0 g (r)V0;(r) - ^(r)V^(r)]. 
The variation of u generated by the external field (73) is 



pq 

= Y\< Ol^Wk > Cu- < v\Ar)\0 > c u ]e- mt 

V 

+ Et< H^)|0 > <£- < 0|J(r)|i/ > cl]e int . (101) 

V 

To proceed further three options are possible. 
6.1 WFM method 

The first way was developed within the WFM approach [6]. It allows one to derive an approxi- 
mate analytical expression for <5u(r, t). The main idea lies in the parametrization of infinitesimal 
displacements £(r, t). Let us recall the main points. By definition 5\ii(r,t) = — — . The 

C6 

displacement & is parametrized [1, 6] by the expansion 

3 3 3 

£i(r,t) = Gi(t) + ^2Gij(t)xj + E Gijk(t)xjX k + E ^jhW^W ^ ( 102 ) 

3=1 j,k=l 3,k,l=l 

which, in principle, is infinite, however one makes the approximation keeping only the first 
terms and neglecting the remainder. For example, in [1] only the two first terms were kept. It 
turned out, that the Gi do not contribute to the final results due to the triplanar symmetry 
of considered nuclei. The coefficients G it j were expressed analytically in terms of the variables 
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1Z 2 i(t) and C n (t). Using the dynamical relations between 72. 21 (t) and £n(t) given by the last 
equation of the set (17), the final formulae for £j(r, t) were found to be 

6 = V2BJ 13 x 3 , 6 = V2BJ 23 x 3 , & = y/2A(J 13Xl + J 23 x 2 ) (103) 

with 

J13 = (£2-1 - ^2l)/2, J23 = i(^2-l + ^2l)/2, 

A= A [ i_ 2 g(l-a)i]/[Qoo(l-|*)], 

5 = ^[1 + 2^(1 - a)5]/[Qoo(l + ^)]. (104) 

By definition the infinitesimal displacements £j are the differentials (£1 = <ix, £2 = c?y, ^3 = dz). 
This fact allows one to construct the differential equations for current fields. For example, for 
the current field in the plane y = we have 

dx B z B 

— = — — — > xdx -zdz = 0. (105) 

dz Ax A v ' 

Integrating this equation we find 

x 2 + az 2 = const = c — > 1 — = 1, 

c c/a 

where o = —B/A. Depending on the sign of a this curve will be either an ellipse or a hyperbola. 
It was shown in [1] that the curve is an ellipse for the scissors mode and it is a hyperbola for 
IVGQR (see Figs 1,2). It was shown also, that the real motion of the scissors mode is a mixture 
of rotational and irrotational behaviour. To get a quantitative measure for the contribution 

— * 

of each kind of motion, it is sufficient to write the displacement £ as the superposition of a 
rotational component with the coefficient a and an irrotational one with the coefficient b [21]: 

£= a ?x x f + bV(yz) = a(0, -z, y) + 6(0, z, y). 

Comparing the components £ y — (b — a)z, (, z = (b + a)y with £ 2 , £3 in (103) we find 

b- a = V2.J 23 B, b + a = V2.J 23 A — ► a = r](l + a), b = r](l-a), 

where n = J 23 A/\^2. So, for the scissors mode in the small 5 limit we have 
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Figure 1: Schematic picture of isovector displacements for the scissors mode. Thin ellipses are 
the lines of currents. The thick oval is the initial position of the nucleus' surface (common 
for protons and neutrons). The dashed oval is the final position of the protons' (or neutrons') 
surface as a result of an infinitesimal displacements shown by the arrows. 
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a = 2r,{\-\s), b=^5, b/a ~ ^(1 + ^5) » (106) 

i.e. the current of the scissors mode is dominated by rotational motion. The contributions of 
the two kinds of motion to the IVGQR are 

a = \v8, 6 = 277(1 - a/b ~ ±<J(1 + ^) » ^, (107) 

i.e. the current of the IVGQR is dominated by irrotational motion. 

Transitions currents are calculated in WFM analogously to transition probabilities. The 
pole structure of the right hand side of equation (101) tells us, that the transition current can 
be calculated by means of an expression similar to (22): 



< 0\Ji(r)\u >= h lim (ft - 0,)p e<? (r)&(r, t) exp(iftt)/ < U \W\0 > . (108) 



For the & from above we obtain (using formulae (77) and (78)) 

71 



A - 
<0|J 3 (r)|z/> = -ifi I/ p c «(r)-=[< 0|-R 2 -i - ^21 k > &i + i < 0|-R 2 -i + #2i|^ > z 2 ], 



<0|J 2 (r)|z/> = 0,^(1-)-/= < 014-1 + 4ik> a*, 

<0|Ji(r)|i/> = -iO,p e<? (rW < 0|i? 2 -i - #21^ > a*- (109) 

v2 

As it is seen, transition currents are proportional to transition probabilities. 

If necessary, one can find the next term of the series (102). To calculate the respective 
coefficients Gijkiit) in the WFM method one is obliged to derive (and solve) the set of dynamical 
equations for higher (fourth) order moments of the Wigner function. Examples of similar 
calculations for third rank tensors can be found in [33]. 

6.2 RPA method 

The procedure of constructing the flow distributions in RPA is more complicated. It is necessary 
at first to calculate transition currents. Having solutions (37) for X^ i: Y£ { one can do it with 
the help of formula (50): 

jim Qim jmi Q-mi 1 



< 0\J(T)\U >= Y,UirnX mi +jmiV mi ) = £ { 
mi mi 

jim Qim jmi Q n 




E v — e E u + eo 



+ E 

mi(AN=2) 
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c mi ~ c rm 

jim Qim jmi Q m 



E v — e 2 E v + e 2 



(110) 



Figure 2: Schematic picture of isovector displacements for the high-lying mode (IVGQR). The 
lines of currents are shown by thin lines (hyperbolae). The thick oval is the initial position of 
the nucleus' surface (common for protons and neutrons). The dashed oval is the final position 
of the protons' (or neutrons') surface as a result of an infinitesimal displacements shown by the 
arrows. 
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The operator Q has a finite number of particle hole matrix elements Q mi , so, in principle, the 
sums in (110) can be calculated exactly. The same is true for the coefficients c„ (78). Therefore, 
in accordance with (101) one could hope to find the exact RPA result for the velocity distribution 
8u(r, t). Unfortunately, because of the pole structure of coefficients c u (Q), it can be done for 
any fl except the required frequency Q u corresponding to the considered mode (resonance). Of 
course, it is clear that in the case of Q close enough to fl v the main contribution into 5u comes 
from the single matrix element < 0| J(r)\u >. That is why, to get an idea about the distribution 
of currents in the RPA eigenstate \v > it is sufficient to know the transition matrix element 
< 0|J(r)|z/ >. However, even in this simple model one can not find a compact analytical 
expression for sums in (110) - the field of velocities can be constructed only numerically. 

As we have already seen it is much more convenient to deal with lines of currents. The 
differential equation for them can be derived with the help of formula (101). With a e~ im time 
dependence we rewrite it in the more convenient form 

-p e «(r)^(r) = ]T[< 0|J(r)|<r > c a - < a\J(r)\0 > c a \. 

<J 

and define the ratio 

6(r) = E CT [< 0|Ji(r)|o- > eg- < o-| Ji(r)|0 > cj 
&(r) E CT [< 0|J 3 (r)|cr > c a - < a\J 3 (r)\0 > c a \ 

Remembering the definition of ^ and c a , multiplying the numerator and the denominator of the 

right hand side by (Q — Q v ) and taking the limit Q — > Q u we arrive to the differential equation 

dx _ <0\J 1 (r)\u> 
dz <0\j 3 (r)\v> 

which determines the lines of currents in the plane y = const for the resonance state \v >. 
6.3 Green's function method 

The distribution function being known, one can calculate the distribution of nuclear currents 
j T (r,t) = mp T (r, t)u T (r, t). There are no any currents in the equilibrium state, so we have 

f x (v,t) = mpT(T,t)6uZ(T,t) = J*?*p x K(T,j>,t) 
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-~[2rne F - nf^ix 2 + y') - m 2 ^Vp[C;(tK + CI(t)w-] 

= -^7^(r)[c;(tv + + c:(tv_], (ii2) 



where the following notation is introduced 

t 

C T ± (t)= J dt'Zl(t') cos u;±(t' -t). 

— oo 

Deriving (112) we used the approximation (11) which means, in particular, that u>f = cuf and 
Po = Po = Po/2. Another component of the flow is 

= / (^3^A T ( r 'P^) = -^^WI^^K " ^V-l, (us) 

With the help of the isovector counterpart of formulae (91) the functions C±(t) can be written 
via dynamical variables 

cz(t) = [L[(t) + ^r y (t)}/(, cut) = -m(t) + —r y (t)}/c 

2n 3 ei 2 

with C = 7 — rr^ — /3mu x uj z /GK,o and the required combinations are 

3mu^ (2nri) 6 

c T + (t)u + + cr_(t)u>- = -2u z [L\(t) - r y (t)]/c, 

Cl(t)u + - C T _(t)u_ = -2u x [L\(t) + F y (t)}/( (114) 

We are interested in isovector flows j x — fx — fx and j z — ft — Jz ■ With the help of the first 
and last equations of (17) we find 

C T + {t)uj + + CT(t)u- = -2u z iQm[l + 2— (1 - a)5}Q 1 /(, 
C T + {t)u + - C T _(t)u_ = -2u x inm[l - 2— (1 - a)5}Q l /C- (115) 
As a result, we have the explicit expressions for the currents 

Ur,t) = ^[l-2^ 2 (l-a)5}Q 1 (t)p (v)^(x + y), 

Ur,t) = ^[1 + 2^(1 -a)8}Q 1 (t)p (v)^z. (116) 
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Following the recipe of section 6.1 we can derive the differential equation for lines of currents, 

for example, in the plane y = 0: 

dx_ _ jx_ dx _ z Jl 1 + 2^(1 - a) 5 _ zB 

dz~] z ~^ dz ~ xuo 2 x 1 - 2g(l - a)5 ~ x A 

with A and B defined by (104). Obviously, this expression coincides exactly with formula 

(105). It is necessary to emphasize the principal point: the result (117) is obtained from the 

GF method in a direct way, whereas deriving formula (105) we made the strong approximation 

about truncating the expansion (102) which parametrizes the displacements. The agreement 

of both expressions tells us about the internal consistency of the various approaches to obtain 

the gross structure of the flow patterns. 

6.4 Summary of flow calculations 

In conclusion in full RPA one must calculate the currents numerically leading to fine details 
(shell effects) whereas in WFM and GF treatments one obtains their gross structure with ana- 
lytical formulas. The latter feature is quite important in order to understand the real character 
of the motion under study since current patterns produced numerically from complicated for- 
mulas with a lot of summations like in (110) can hardly be interpreted physically. A good 
example is the interplay of the scissors mode and the isovector giant quadrupole resonance. 
Looking only at the flow patterns (see Figs. 1, 2) one would not be able to tell that the former 
is mostly rotational with a small amount of an irrotational component and the other way round 
for the latter, as this can be seen from eqs. (106, 107). 

7 Conclusion 

In this paper we made an exhaustive comparison of different methods to treat collective excita- 
tions in nuclei, like the scissors mode, isovector and isoscalar giant quadrupole resonances. This 
comparison was exemplified on the H.O. plus separable quadrupole-quadrupole force model but 
it has more general character. 

We investigated WFM, RPA and Green's Function (GF) methods. Under certain circum- 
stances all three methods give essentially the same results. For example all methods give in 
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our model the same analytical expressions for energies and transition probabilities for all the 
excitations considered. It turned out that WFM and GF methods are very close to one an- 
other. Contrary to RPA, both work in phase space and incorporate semiclassical aspects, with 
no need to introduce a single particle basis. Finally both methods yield identical sets of dynam- 
ical equations for the moments. However, in the case of realistic forces the GF method loses 
its simplicity and the more complicated pseudo particle method has to be applied, whereas 
the WFM method continues to be a convenient and powerful tool for the description of the 
collective motions as it was demonstrated in ref. [6, 25, 26] employing Skyrme forces. 

To show the analytical equivalence between WFM and RPA methods one needs to introduce 
the dynamical equations for the transition matrix elements. They can be derived either from 
the RPA equations for the amplitudes X kq , Y kq or from the WFM dynamical equations for the 
moments. This proves the identity of eigenvalues in both methods under the condition that 
a complete basis is used in both cases. However, both methods behave differently when the 
dimension of the space is reduced. Actually WFM is designed to use only rather few moments of 
low rank. The restricted number of eigenvalues approximate the collective states in an optimal 
way, representing e.g. their centroid positions, as this was shown in [6, 25, 26, 27]. In this sense 
WFM has similarity with the sum rule approach [4] which works, however, only in the cases 
when practically all strength is exhausted by one state, whereas WFM method works also in 
situations when the strength is distributed among several excitations. On the contrary in RPA 
one needs in general a rather large space to correctly account for the collectivity of e.g. the 
giant resonances. At the same time a certain fine structure of the resonances is also obtained. 
Concerning the spectrum both methods are thus complementary. The situation is different for 
the currents and flow patterns. Since RPA is a fully quantal approach, the current lines can 
even in our simplified model be calculated only numerically. They show a complicated pattern 
due to the shell effects. Without further efforts one will not be able to analyze the nature of 
the flows. A good example is given in our model with the low and high lying scissors mode 
(the latter being the IVGQR). Due to the analytic formulas found with WFM and GF methods 
which naturally lead to smooth current distributions free of rapidly fluctuating behavior from 
shell effects, we were able to show that the low lying scissors mode is mostly rotational with a 
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slight admixture of an irrotational component and the other way round for the IVGQR. 

In addition to our earlier work we investigated in detail the so called synthetic scissors mode 
which is based on the picture of two counter rotating proton and neutron mass distributions. 
Calculating the overlap of this synthetic scissors mode with the real one we could show that the 
squared overlap amounts only to about 60% in the best of all cases. We also showed explicitly 
the orthogonality of the spurious mode to all other "intrinsic" excitations of the model. 

Future work in this direction shall be concerned with the scissors mode in neutron rich 
nuclei and with the consideration of superfluidity. 



Appendix A 

It is known that the deformed harmonic oscillator Hamiltonian can be obtained in a Hartree 
approximation "by making the assumption that the isoscalar part of the QQ force builds the 
one-body container well" [10]. In our case it is obtained quite easily by summing the expressions 
for V p and V n (formula (7)): 

l/(r 1 t)4(P(r,t) + r(r,t)) = ^V + Ko £ (-l)^-^)^)- (A.l) 

1 1 /i=-2 

In the state of equilibrium (i.e., in the absence of an external field) Qi±\ = Q2±2 = 0. Using 
the definition [31] Q 2 o = Qoo^S and the formula g 20 = 2z 2 — x 2 — y 2 we obtain the potential of 
the anisotropic harmonic oscillator 

nr) = ^(* 2 + y 2 )+^V] 

with oscillator frequencies 

uj 2 x = L0 2 y =to 2 (l + a5), u 2 z =u 2 (l-2a5), 

where a = —k , The definition of the deformation parameter 5 must be reproduced by 
3mu 2 

the harmonic oscillator wave functions, which allows one to fix the value of a. We have 

Qw = -{- + ^ + -l Q 20 = 2-(^-^), 

m uj x ojy lo z m oo z oo x 
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where T, x = ^f =1 (n x + and n x is the oscillator quantum number. Using the self-consistency 
condition [31] 

^x^x — ^y^y = ^z^z — ^Wo, 

where S and ujq are defined in the spherical case, we get 

Q20 _ uj 2 x - m 2 2a8 4 
Qoo oo 2 + 2uj 2 l-a5 3 ' 

Solving the last equation with respect to a, we find 

"=3Tm- (A - 2) 

Therefore, the oscillator frequences and the strength constant can be written as 

<f x = <f v =Q\l + \s), ^ z =^(l- 2 -5), ko = ~ (A.3) 

with uj 2 = u 2 /(l + The condition for volume conservation uj x uj y uj z = const = Uq makes u> 
(^-dependent 

2 _ 2 1 + \5 

^ ~ ^(1 + |5) 2 /3(1 _ |^)l/3- 

So the final expressions for oscillator frequences are 



^=^=^^ ■ ^°UrfJ ' (A ' 4) 

It is interesting to compare these expressions with the very popular [31, 4] parametrization 



The volume conservation condition gives 



12 _ U Q 



UJ = 



(i + f^ / ) 2/3 (i — 1^0 1/3 ' 

so the final expressions for oscillator frequences are 

<^-. 2 (^ff > ^(Trff ■ (A ' 5) 

The direct comparison of expressions (A. 4) and (A. 5) allows one to establish the following 
relation between 5 and 5': 

I .6' 



1 + 25' 1-25' 
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One more parametrization of oscillator frequences can be found in the review [21]: 



, .112 112 
2 2 U 2 U 



One has from the volume conservation condition 



^ 2 = c 2 (l-^") 2/3 (l + ^") 1/3 , 
so the final expressions for oscillator frequences are 

2 2 2 A + f<n 1/3 2 2 /l-|^\ 2/3 , A x 

^=«:=^( I r^J . ^=^(tt|^J • (A - 6) 

that coincide exactly with (A. 4), i.e. 5" = 5. 

It is easy to see that equations (A. 4) correspond to the case when the deformed density p(r) 
is obtained from the spherical density po( r ) by the scale transformation [9] 

(x,y,z)^(xe a / 2 ,ye a / 2 ,ze- a ) 

with 

e «= [i±4l] , -, (A.7) 

which conserves the volume and does not destroy the self-consistency, because the density and 
potential are transformed in the same way. 

It is necessary to note that Q 00 also depends on 5 

ft , Xj t ^j?/ f^' / 2 1 . rt X 



^oo-— (— + — + — ) - — N: "^ ll! ~ + ~a> - 



Ww u;,' m u ^ + 1^)1/3(1 -15)2/3' 

where Q Q m = A^R 2 , R = r^A 1 ^. As a result, the final expression for the strength constant 
becomes 

mw 2 ( _ mu 2 _ a 

K °~ 4Q8o + " 4g° e ' 

that coincides with the respective result of [9]. 



Appendix B 

To calculate the sums Q = |Q m i| 2 and Q 2 = V] |Q m i| 2 we employ the sum-rule 

mi(AN=0) mi(AN=2) 

techniques of Suzuki and Rowe [9]. The well known harmonic oscillator relations 



h 



n x +l) t 



y x 
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„ , . / mfiuj x , , , . . N 

PxV'n* = -n/ z {VnxiPn x -l ~ Vn x + llpn x +l) (B.l) 



allow us to write 

ft 



xzip nx tp nz = ■=={y/n x n g tl) nx - 1 tl) nt -i + J (n x + l)(n z + l)ip nx+1 ip nz+1 



+ y (n x + 1)^^+1^,-1 + \/ n ^( n ^ + l)^n„-l^n,+l), 



-Vw0n, = -~ -f={y/n x n z i} nx - 1 i} nz - 1 + ^(n a . + l)(n* + l)V>n ie +i^ 



n z +l 



m z uj x uj z 2,m^/uj x uj z 

-yj (n-x + l)n z ^n B +i^n,-i - \Jn x (n z + 1)^-1^+1)- (B.2) 
These formulae demonstrate in an obvious way that the operators 

Pq = \(zx^ p x p z ) and P 2 = \{zx — ^ p x p z ) 

I m z uj x uj z 2 m z uj x uj z 

contribute only to the excitation of the AiV = and AiV = 2 states, respectively. Following 

[9], we express the zx component of r 2 Y 2 \ = \J~j^Q = — \f^z(x + iy) as 

zx = P + P 2 . 

Hence, we have 

A A 

e o Yl I <0\J2z s x s \mi> | 2 = e ]T| < 0\J2 Po(s)\mi > | 2 

mi(AN=0) s=l mi s=l 

= ^<0|EPo(s),[^i;Po(s)]]|0>, (B.3) 

Z s=l s=l 

where eo = ftf^ — a; z ). The above commutator is easily evaluated for the Hamiltonian with 
the potential (A.l), as 

< 0|EP ( S ), [#,£P (s)]]|0 >= -^0 f <0|Ell ^ |0> - <0|EllX ' |0> ) • (B.4) 



Taking into account the axial symmetry and using the definitions 

A A 

^(Z 2 -X 2 )\0>. Onn = Onr, 

S J. 

we transform this expression to 



Qoo =< 0| £(2:r 2 + z s)\° >> Q20 = 2 < 0| J> 2 - a; 2 )|0 >, Q 20 = qJ^S, 

s=l s=l 



< 0\[^P (s),[H,^Po(s)]]\0 >= —e Q 00 [—^ 8- . (B.5) 
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With the help of the self-consistent expressions for u x , ui z (A. 3) one comes to the following 
result: 

< 01 E \P, t A.W1H0 >= f^l = £QSo (B.6) 

By using the fact that the matrix elements for the zy component of r 2 F 21 are identical to those 
for the zx component, because of axial symmetry, we finally obtain 

By calculating the double commutator for the P 2 operator, we find 

V" I <r- nl Y% 2 v l™; I 2 - 5 r 2 - 5 e2 f 1 + i S \ 1 

e 2 2^ < U 2^*21 ™ > — — T^e 2 - — 2 i 27 ' 

mi(Z^=2) s=i 16vrmw 2 16tt m cu 2 v l-§5/ 



(B.8) 



where e 2 = h(u x + u; z ). 

We need also the sums Qq and <2 2 calculated separately for neutron and proton systems 
with the mean fields V n and V p , respectively. The necessary formulae are easily derivable 
from the already obtained results. There are no reasons to require the fulfillment of the self- 
consistency conditions for neutrons and protons separately, so one has to use formula (B.5). A 
trivial change of notation gives 

z z h /l + ^ p 1 — 2^p\ 
< 0\[Y { Po(s),[H p ,Y { Po(s)]]\0 >= J")' ( R 9) 

e P o E l<0|E^2iM>| 2 = ^^ e ^ (^f^-^-^), (B.10) 



mi(AN=0) s=l 









1+§5P 


1 




UJx 




u? z ) 


1 + |5P 


1 

+ - 





4 E l<0|Er s V 21 H>| 2 = ^A e PgP o (^^ + ^L). (B.H) 

mi(AN=2) 8=1 IDTTm \ UJx UJ Z J 

The nontrivial information is contained in oscillator frequences of the mean fields V p and V 11 
(formula (7)) 

(cup) 2 = uj 2 [1 - -^(kQ» + RQ» )}, K) 2 = u 2 [l + — 2 {kQI, + RQ^)\, 



(cu») 2 = u 2 [l - -^(kQSo + RQ p 20 )}, K) 2 = uj 2 [1 + — 2 {kQI q + «<&)]. (B.12) 



48 



The above-written formulae can also be used to calculate the analogous sums for operators 
containing various combinations of momenta and coordinates, for example, components of an 
angular momentum, tensor products (rp) 2 i and (j3 2 ) 2 i. By definition I 1 = yp z — zp y , I 2 = 
zp x — xp z . In accordance with (B.l), we have 



XPz^n x lpn z = -^J — (V n x n z^n x -l^n z -l ~ \J ijl x + \){n z + 1)^+1^ 



n z +l 



+ \j (n x + 1)^^+1^-1 - \Jn x (n z + l)Vn x -iVv+i)- 



(B.13) 



Therefore, 



hipn x tpn z = J— - x \ — ){^/n x n z i) nx -ii) nz -i - \J{n x + l)(n z + l)V>n x +iVv+0 



UJ X 



+i oUl + + 1)^^+1^,-1 - \jn x {n z + l^-iVv+i)- (B- 14 ) 

Having formulae (B.2) and (B.14), one derives the following expressions for matrix elements 
coupling the ground state with AiV = 2 and AN = excitations: 



< n x + 1, n z + 1|/2|0 >= i 



: h(uj 2 x -uj 2 z ) (n x + l)(n z + l) 



2 uj x + uj z 



Ld x LU z 



< n x + l,n z — 1 1 ^2 1 >= i 



< n x + 1, n z + l|a;^|0 >= 



2 u; T 



"V 



h (n x + l)(n z + l) 



2m 



x^z 



■ 1 ft / (n^ + l)n z 
< n x + l,n z — l\xz\0 >= \' 



2m 



LU X LU Z 



(B.15) 



It is easy to see that 

< n x + 1, n z + 1|/2|0 >= im- 



"1 - u 2 z ) 
uj x + u z 



< n x + 1, n z + l|x2;|0 >, 



< n x + l,n z - l\i 2 \0 >= im— — < n x + l,n z - l\xz\0 > . 

uj x — u z 

Due to the degeneracy of the model all particle hole excitations with AiV = 2 have the same 
energy e 2 and all particle hole excitations with AiV = have the energy eo- This fact allows 
one to join the last two formulae into one general expression 



< ph\l2\0 >= ihm- 



€ph 



< ph\xz\0 > . 
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Taking into account the axial symmetry we have an analogous formula for 1\. 

(lu 2 — OJ 2 ) 

< ph\h\0 >= -ihm— — < ph\yz\0 > . 

The magnetic transition operator (29) is proportional to the angular momentum: / 1±1 

ie I 3 - 

— (1 2 =F Hi) Therefore, we can write 



4mc V 27r 

< p/»|/i±i|0 >= e M"« ~ "*) < ph\r 2 Y 2±1 \0 > . (B.16) 
2c V 5 £ P h 

Similar calculations for the tensor product (rp) 2 i = —\[zp x + xp z + i(zp y + yp z )\ lead to the 
following relation: 

Tfl / 27T Tfl 

< ph\(rp) 2 i\0 >= —e ph < ph\r 2 Y 2±1 \0 >= i^ P h < P%2il° > • ( B - 17 ) 
Two kinds of particle hole matrix elements are obtained from the second formula of (B.2): 



< n x + l,n z + l\pxPz\Q >= —Timuj x u) z 

<n x + l,n z — l\p x p z \Q >= hmuj x uj z \ 
Simple comparison with (B.15) shows that 



'(ra* + l)K + l) 



2u; T 2u;, 



! (n x + l)n z 
2u x 2uj z 



< n x + 1, n z + llPsPziO >= —m 2 uj x uj z < n x + 1, n z + l|rr^|0 >, 

< n x + 1, n z — l\p x p z \0 >= m 2 uj x uj z < n x + 1, n z — l\xz\0 > . 
With the help of the obvious relations 

2u) x u) z = u) 2 x + uj 2 z - e 2 /h 2 , -2u) x u) z = u 2 x + u 2 z - e\jU 2 

these two formulae can be joined into one expression 



< ph\p x p z \0 >= — (oj 2 x + uj 2 - e 2 ph /h 2 ) < ph\xz\0 > . 
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By definition p\ x = —p z {p x + ip y ) and f\ x = —z{x + iy), hence, 



m 



2 



< Ph\p 2 21 \0 >= — K + < ~ 4h/ h ) < P h \ r 2i\° > ■ (B.18) 
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